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Detsmina^-ion of the correlation phase between a signal 
and a replica sequence 

FIELD OF THE INVENTION 

The invention relates to a method for determining the 
correlation phase between a signal received at a receiver 
and an available replica sequence by using a matched 
filter checking various correlation phases. The invention 
relates equally to a receiver, to an electronic device 
and a communication system comprising such a receiver, 
and to a device communicating with such a receiver. 

BACKGROUND OF THE INVENTION 



The correlation phase between a signal received at a 
receiver and an available replica sequence has to be 
determined for example for CDMA (Code Division Multiple 
Access) spread spectrum receivers. 

For a spread ..spectrum communication in its basic form, a 
data sequence is used by a transmitting unit to modulate 
a sinusoidal carrier and then the bandwidth of the 
resulting signal is spread to a much larger value. For 
spreading the bandwidth, the single- frequency carrier can 
be multiplied for example by a high-rate binary pseudo- 
random noise (PRN) code sequence comprising values of -1 
and 1, which code sequence is known to a receiver. Thus, 
the signal that is transmitted includes a data component, 
a PRN component, and a sinusoidal carrier component. A 
PRN code period comprises typically 1023 chips, the term 
chips being used to designate the bits of the code 
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bits of the data sequence. 
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a parallel processing for all possible correlation 
phases. Thus, the DFT based correlators enable a 
reduction of the computational burden, in particular if 
the DPTs are realized as Fast Fourier Transforms (FFT) . 

Nevertheless, the first type of correlators is still 
widely used. One reason for this is the simplicity of the 
implementation. While the FFT might be suited to reduce 
the complexity of the correlator, the use of FFT requires 
complicated algorithms. Moreover, the FFT approximates 
calculations due to the use of transform coefficients. 
Another reason for further performing calculations in the 
time domain is that while in principle, the phase of the 
received signal relative to the available replica 
sequence can have any possible value, in some situations, 
the range of the possible correlation phases can be 
reduced based on some a-priori knowledge. Currently, 
however, a limited search of correlation phases can only 
be realized with correlators performing a correlation in 
the time domain. Known DFT based methods inherently 
perform the search of all possible correlation phase in 
parallel- Therefore, their usage is not feasible in 
situations in which the search is to be carried out only 
over a limited number of all possible correlation phases. 

SUMMARY OF THE INVENTION 

It is an object of the invention to reduce the 
computational load of correlators performing a 
correlation in the time domain. More specifically, it is 
an object of the invention to reduce the computational 
load of matched filters which are employed for 
determining the correlation phase between a received 
signal and an available replica sequence. 
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In a method according to the invention, a matched filter 
multiplies samples of the received signal with samples of 
the replica sequence and sums the resulting products to 
obtain a correlation value for a specific correlation 
phase. The multiplication performed by the matched filter 
can be for instance an elementwise multiplication. The 
samples of the received signal and the samples of the 
available replica sequence are shifted relative to each 
other for each correlation phase which is to be checked. 
It is proposed that results obtained in the calculations 
for one correlation phase are used by the matched filter 
also for calculations for a subsequent correlation phase. 

A conventional matched filter, in contrast, performs the 
calculations for each correlation phase independently. 

Moreover, a receiver, an electronic device comprising a 
receiver and some other device are proposed, either 
comprising means for carrying out the steps of the 
proposed method. The electronic device can be for 
instance a mobile terminal which is capable of 
communicating with a communication network. In case the 
processing is performed in another unit than the 
receiver, the required information about the received 
signals is forwarded by the receiver to this unit. The 
proposed other device can be for instance a network 
element of a network. The object is also reached with a 
system comprising a receiver and a device, in which 
system either the receiver or the device comprises means 
for carrying out the steps of the proposed method. In 
case the receiver performs the processing, the device may 
provide assistance data to the receiver. 
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The invention proceeds from the consideration that 
usually, a replica sequence will comprise several sub- 
sequences of at least twice the same value in a row. When 
the received samples and the replica samples are shifted 
by one position relative to each other after a 
correlation value for one correlation phase has been 
determined, thus partly the same multiplication results 
are required again when determining the correlation value 
of the next correlation phase. This occurrence can be 
exploited for reducing the number of multiplications 
and/or summings for a respective next correlation phase, 
since the used replica sequence is known. 

It is an advantage of the invention that the number of 
operations can be reduced, and thus the employed 
processing means, either hardware or software, may be 
slower compared to known methods. The frequency 
requirement is reduced for software, which allows to 
employ software correlators with a lower frequency. Thus, 
the invention facilitates the use of software 
correlators. If the operations are carried out slowly, 
also the power consumption can be reduced. 

Alternatively, the reduced number of operations can be 
used for carrying out the processing faster than before. 
If the operations are carried out fast, the acquisition 
of signals and following operations like a positioning 
can be accelerated. 

Due to the reduced number of operations, the matched 
filter is also particularly easy to implement and more 
competitive with DFT based correlators. 
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Preferred embodiments of the invention become apparent 
from the dependent claims. 

The amount of operations can be reduced with the 
invention for any kind of signal which has to be 
correlated with an available replica sequence. 
The received signal can be for instance a code modulated 
signal and the replica sequence a replica code sequence. 

The invention is of particular advantage, when the 
received signal comprises a binary sequences, for 
instance a BPSK (binary phase shift keying) binary code 
with which the signal is modulated. A BPSK binary code is 
employed e.g. in GPS. A binary sequence is composed of 
the values +1 and -1. In case the signal comprises a 
binary sequence, only about half of the elements will 
change their values from one replica sample to the next. 
The required amount of operations for the correlation can 
therefore be halved with the invention. As a result, the 
hardware or software could be half as fast to obtain a 
correlation phase with the same speed as with a known 
matched filter, or the acquisition of a correlation phase 
could be twice as fast compared to a known matched filter 
when using the same hardware or software. Still, the 
invention is also of advantage when used in other 
systems, in which the two input signals to the matched 
filter may be composed of any number of any discrete 
values . 

The multiplications in the correlation calculations can 
be performed in different ways. For example, in case the 
received signal and the available replica sequence 
comprise a binary sequence consisting of the values +1 
and -1, the multiplications may be realized e.g. simply 
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with a sign change for each required multiplication by 
-1. 

The results of correlation calculations for a preceding 
correlation phase can be made use of in several ways when 
determining the correlation value for the current 
correlation phase. The most suitable way may also depend 
on the kind of sequence comprised in the received signal, 
e.g. on the employed code. For example, the entire 
correlation value resulting for a preceding correlation 
phase can be used as basis, which is then corrected for 
the current correlation phase. In one alternative, only 
those multiplication results of a preceding correlation 
phase which are known to stay the same are re -used for 
the current correlation phase. 

The invention can be used in both, acquisition and 
tracking schemes. 

In both cases, the invention may be used for determining 
the correlation phase and the frequency of a remaining 
complex sinusoidal modulation. The correlation phase is 
determined according to the peaks of a cross -correlation 
function, and the correlation is calculated e.g. at 
initial code wipe-off stages. The processing for weak 
signals requires additional coherent and non- coherent 
integrations. The invention can therefore also be used as 
a building block for other methods implementing different 
scenarios of coherent and/or non-coherent processing for 
possible multiple frequency candidates. 

As indicated above, the invention can be implemented in 
hardware or in software. In case the invention is 
employed as part of acquisition and tracking algorithms, 
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the implementation corresponds advantageously to the 
implementation of these algorithms. 

The invention can be used in any system which uses 
matched filters. It can be employed in particular, though 
not exclusively, in a communication system using CDMA 
spread spectrum receivers, for instance in a positioning 
system like GPS or Galileo. Alternatively, it could be 
employed for example in radar technology. 

BRIEF DESCRIPTION OP THE FIGURES 

Other objects and features of the present invention will 
become apparent from the following detailed description 
considered in conjunction with the accompanying drawings, 
wherein: 

Fig. 1 illustrates an embodiment of the method according 

to the invention; 
Fig. 2a shows a sequence of samples of an received code 

modulated signal; 
Fig. 2b shows samples of a binary replica code sec[uence; 

and 

Fig. 2c indicates sign changes in the replica code 
sequence of figure lb. 

DETAILED DESCRIPTION OF THE INVENTION 

Figure 1 illustrates an exemplary embodiment of the 
method according to the invention implemented in a GPS 
receiver. 

The GPS receiver comprises a receiving unit for receiving 
signals from different GPS satellites which are modulated 
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with different C/A-codes. The C/A codes employed by the 
satellites for modulating the signals are BPSK binary PRN 
codes consisting of only +1 and -1. The GPS receiver 
moreover comprises a tracking unit with a correlator. The 
correlator employs a matched filter for acquiring and 
tracking received satellite signals. The depicted method 
is implemented more specifically in this matched filter. 
The receiver further stores a replica code sequence for 
each of the satellites, which are also called space 
vehicles SV. Also the positions of sign changes in the 
stored replica code sequences are known to the receiver. 

For determining the correlation between a signal received 
from one of the satellites SV and the corresponding 
replica code sequence, first this replica code sequence 
is retrieved. The number of samples of the replica code 
sequence (N) depends on the length of the code employed 
by the respective satellite for modulating the signal. It 
also depends on the sampling rate of the incoming signal. 
E.g., in case of 1 sample per chip sampling, the replica 
code sequence length is W=1023. The samples of the 
retrieved replica code sequence are denoted as: 

A counter value i counting the code phases that are to be 
checked is set to i=0. 

Further, a first set of consecutive samples of a received 
code modulated signal is stored in a memory of the 
receiver. The received and stored samples are denoted as: 

^N-2f • " • / "^3' ^1' ^0 
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These samples are multiplies elementwise with samples of 
the retrieved replica code sequence. 

The resulting products are then summed, in order to 
obtain a correlation value Co for the first code phase. 

In the following, a loop is entered for determining the 
correlation values also for subsequent code phases. For 
each code phase that is to be checked, the counter value 
i is incremented by 1 at the beginning of the loop. 

Further, the memory of the receiver is updated for each 
code phase by removing the respective oldest sample and 
by introducing a new sample. For the second code phase, 
for example, a new sample of the incoming signal is 
added to the set of consecutive samples, while the oldest 
sample xo is removed from the memory. The stored set of N 
samples for the second code phase is thus given by: 

The replica code sequence, in contrast, stays the same 
for all code phases, therefore a relative shift is 
obtained between the received samples and the replica 
samples . 

Based on the respective new set of code samples, a 
correlation value Ci is also calculated for the second 
code phase and later on for the subsequent code phases. 
In contrast to a conventional correlation, however, 
multiplications are not carried out for all samples of 
the respective set of received code samples, as will be 
explained in the following. 
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When comparing the first set of stored input samples with 
the second set of stored input samples, it can be seen 

that the samples jca^_i,x;v^_2»— »^3»-^2>-^l occur in both sets-, 
even though they are aligned differently with the replica 
samples. In case a shifted received sample is aligned 
with a replica sample of the same value as before, the 
multiplication for this sample will also have the same 
result as before, which can be made use of. 

The conventional correlator output, i.e. the correlation 
value Ci for a specific code phase i, can be written as: 

N-1 

First this sum can be split, taking out the respective 
oldest sample x^i 

N-1 N-1 

Then, the remaining sum can be considered to. be composed 
of two subsums : 

The first subsum contains the multiplications of input 
samples with those replica samples which have the same 
sign as the respective preceding replica sample. The set 
of the values of j identifying these replica samples is 
denoted as J+, where 



- 13 - 



J € J, if = 1 & rj., = l) OR (rj = -1 & r,., = -l) . 

The second subsum contains the multiplications of input 
samples with those replica samples which have another 
sign than the respective preceding replica sample. The 
set of the values of j identifying these replica samples 
is denoted as J., where 

D € J", if (r^ = 1 & r^.^ = -l) OR (r^ = -1 & r^.^ = l). 

For illustration, figure 2a shows a row of eleven 
rectangles 21 representing consecutive samples of a part 
of a stored set of input samples without indication of 
any specific values. Figure 2b shows a row of eleven 
rectangles 22 representing samples of a part of an 
exemplary replica code sequence having from left to right 
the values {+1, -1, +1, -1, -1, +1, +i, +i, +i, -i, -i}. 
The indices of the input samples x and the replica 
samples r increases from right to right. While the 
replica samples remain unchanged for all code phases, the 
input samples are shifted for each new code phase by one 
sample to the right, such that the oldest sample is 
removed from the set of samples and a new sample is 
introduced to the set of samples at the left hand side. 

Figure 2c points out by means of a further row of eleven 
rectangles 2 3 those replica samples of figure 2b which 
provide the same sign for an associated input sample of 
figure 2a as will be provided to this input sample after 
a shift, i.e. those replica samples rj which are identical 
to the respective preceding replica sample rj.i. These 
replica samples are indicated with an associated white 
rectangle, while the replica samples providing another 
sign than the preceding replica sample are indicated with 
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an associated black rectangle. Thus, the 4^^, 6^^ , 7^^, 8^^ 
and 10^^ rectangle from the left is white, and the 1^*^, 
2^*^, 3^^, 5*^^ and 9^^ rectangle from the left is black. 
Since only part of the entire code is depicted, the 11^^ 
replica sample from the left cannot be told from the 
figure to repeat or invert the preceding sample, and thus 
a hatched rectangle is associated to it. 

The elementwise multiplication with' the replica samples 
indicated with a white rectangle will result in the same 
product before a shift as the elementwise multiplication 
with the preceding replica samples after a shift. For 
example, if the 6^^ input sample from the left is assumed 
to have a value "a", the elementwise multiplication will 
result in (+l)*a, since the value of the 6^^ replica 
sample is +1. After the shift by one sample to the right 
for the next code phase, "a" will constitute the value of 
the 7^^ sample from the left in figure la. Since the 7^^ 
replica sample has equally a value of +1, the elementwise 
multiplication will result again in (+l)*a. 

The indices of those replica samples to which a white 
rectangle is associated are part of the set tJ+, while the 
indices of those replica samples to which a black 
rectangle is associated are part of the set J.. 

With the above definitions, the correlation value Ci+i for 
a respective next code phase i+1 can be written as: 



or alternatively as: 
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In the example presented in figures 2a to 2c, the sets J* 
and J. are approximately of the same size. This means that 
with both equations, instead of performing N additions at 
each stage, only N/2+2 additions are needed, plus some 
shifts for realizing the multiplications by 2. These 
shifts could be incorporated into the addition operation 
without increasing the cost of operation. 

Sometimes, however, the distribution of the indices of 
the replica samples to the sets J+ and J. may vary between 
different replica code sequences. 

In the embodiment of the method illustrated in figure 1, 
it is therefore first checked, which one of the sets ^7+ 
and J- has a larger size. For each replica, the equation 
(1) or (2) is selected which causes the minimum number of 
operations depending on the number of samples in the sets 
J+ and J"-. That is, in case the set iJ+ comprises less 
samples, equation (1) is used, while in case the set J. 
comprises less samples, equation (2) is used for 
determining the correlation for the next code phase. 
Thus, there will be different equations used for 
different codes. 

It is to be noted that the actual comparison between the 
sizes of sets J"+ and J*, does not have to be carried out in 
real time during the matched filter operation. Since the 
replica code sequences for all satellites are known and 
the sign-changes within the replica code sequences as 
well, the sizes of the two sets can easily be determined 
and compared beforehand and the result be stored for use 
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by the matched filter. Then, it will be known at the 
matched filter without any computational requirements 
that for one satellite the first option should be used, 
while for another satellite, the secuiid option should be 
used. With this approach, at least a reduction in 
complexity by half is achieved. 

Alternatively, one of the equations (1) or (2) could be 
detezTTiined to be used for all replica code sequences. In 
this case, a reduction in complexity by half will be 
achieved on an average. 

The described loop will be continued for a respective 
subsequent code phase, until a correlation value C± has 
been determined for all N code phases that are to be 
checked. 

The correct code phase can then be obtained by 
determining the maximum correlation value out of all 
determined correlation values. 

The presented embodiment can easily be generalized for 
replica sequences composed of samples with values from 
any discrete set, preferably with a small number of 
elements. For binary replica sequences, this set contains 
only +1 and -1, as mentioned above, and has thus just two 
elements. 

It is to be noted that the described embodiment 
constitutes only one of a variety of possible embodiments 
of the invention. 
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Claims 

1. Method for determining the correlation phase between 
a signal received at a receiver and an available 
replica sequence by using a matched filter checking 
various correlation phases, said matched filter 
multiplying samples (21) of said received signal with 
samples (22) of said replica sequence and summing the 
resulting products to obtain a correlation value for 
a specific correlation phase, which samples (21) of 
said received signal and which samples (22) of said 
available replica sequence are shifted relative to 
each other for each correlation phase which is to be 
checked, wherein results obtained in the calculations 
for one correlation phase are used by said matched 
filter for calculations for a subsequent correlation 
phase . 

2. Method according to claim 1, wherein said matched 
filter multiplies said samples (21) of said received 
signal elementwise with samples (22) of said replica 
sequence . 

3. Method according to claim 1 or 2, wherein said 
received signal comprises a binary secjuence. 

4. Method according to claim 3, wherein the possible 
values of said binairy sequence are +1 and -1. 

5. Method according to claim 4, wherein Ci constitutes a 
determined correlation value for an i*^*^ checked 
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correlation phase, wherein N is the length of said 
binary sequence, wherein Xj+i constitutes a j^^ samples 
of said received signal for said i^^ correlation 
phase, wherein the aaruples (22) of said replica 
sequence are not shifted for the different 
correlation phases which are to be checked, wherein rj 
constitutes a j^^ sample of said replica sequence, 
wherein a set J+ comprises the indices j for which 
= 1 & r,., = l) OR (r, = -1 & r,,,* = -l), 

and wherein a correlation value Ci+i for the (i+1)^^ 
correlation phase is calculated as: 

Method according to claim 4, wherein Ci constitutes 
the determined correlation value for an i*^^ checked 
correlation phase i, wherein N is the length of said 
binary sequence, wherein Xj+i constitutes a j*^^ samples 
of said received signal for said i^^ correlation 
phase, wherein the samples (22) of said replica 
sequence are not shifted for the different 
correlation phases which are to be checked, wherein rj 
constitutes a j^^ samples of a said replica sequence, 
wherein a set J. comprises the indices j for which 
(r,. = 1 & r.., = -l) OR {r^ = -1 & r,_, = l), 

and wherein a correlation value Ci+i for the (i+1)^^ 
correlation phase is calculated as: 

Method according to claim 4, wherein Ci constitutes a 
determined correlation value for an i^*^ checked 
correlation phase, wherein N is the length of said 
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binary sequence, wherein Xj+i constitutes a j^*^ samples 
of said received signal for said i^^ correlation 
phase, wherein the samples (22) of said replica 
sequence are not shifted for the different 
correlation phases which are to be checked, wherein rj 
constitutes a j^^ sample of said replica sequence, 
wherein a set J+ comprises the indices j for which 
k = 1 & r,., - l) OR (r, = -1 & r^.., = -l), 
wherein a set J- comprises the indices j for which 
(r^. = 1 & r^^, = -l) OR (r,. = -1 & r^.., = l), 

and wherein a correlation value Ci+i for the (i+1)^*^ 
correlation phase is calculated as : 

if the size of said set J. is larger than the size of 
said set J^, and as: 

if the size of said set J+ is larger than the size of 
said set J. . 



8. Method according to one of the preceding claims, 
further comprising a subsequent coherent and/or 
noncoherent processing for handling signals of low 
strength. 

9. Method according to one of the preceding claims, 
wherein said received signal is a code modulated 
signal, and wherein said replica sequence is a 
replica code sequence. 

10. Method according to claim 9, wherein said code 
modulation of said received code modulated signal is 
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a Code Division Multiple Access (CDMA) spread 
spectrum modulation. 

11. Use nf a rn^thod ciccording to one of the preceding 
claims in a process for acquisition and/or tracking 
of signals received at a receiver. 

12. Receiver comprising 

receiving means for receiving signals; and 
- processing means for carrying out the method 
according to one of claims 1 to 10. 

13. Receiver according to claim 12, which receiver is a 
receiver of a positioning system. 

14. Electronic device comprising a receiver according to 
claim 12 or claim 13. 

15. Electronic device according to claim 14, wherein said 
electronic device is a mobile terminal capable of 
communicating with a communication network. 

16. Device comprising 

means for receiving from a receiver information on 
signals received by said receiver; and 
processing means for carrying out the method 
according to one claims 1 to 10. 



17. Device according to claim 16, which device is a 
network element of a network. 



System comprising 

a receiver comprising means for receiving signals, 
and means for providing information on received 
signals : and 

a device according to claim 16 or 17. 

System comprising 

a receiver according to claim 12 or 13; and 
a device for providing assistance data to said 
receiver. 

System according to claim 19, wherein said device is 
a network element of a network. 

System according to one of claims 18 to 20, wherein 
said system is a positioning system. 
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Abstract 

The invention relates to a method for determining the 
correlation phase between a signal received at a receiver 
and a replica sequence. A matched filter multiplies 
samples 21 of the received signal with samples 22 of the 
replica and sums the resulting products to obtain a 
correlation value for a specific correlation phase. The 
samples of the received signal and the replica are 
shifted relative to each other for each correlation phase 
that is to be checked. In order to reduce the 
computational load, it is proposed that results obtained 
in the correlation calculations for one correlation phase 
are used by the matched filter also for calculations for 
a subsequent correlation phase. The invention relates 
equally to a corresponding receiver, to an electronic 
device comprising such a receiver, to a device 
cooperating with such a receiver and to a corresponding 
system. 
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